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Abstrsct 


The  asymptotic  distribution  of  ths  log-likellhood  ratio  Is 
shown  to  provide  a  method  of  determining  approximate  confidence  limits 
for  any  coherent  system  when  each  component  has  an  exponential  life 
with  unknown  failure  rate  and  component  performance  data  are  provided 
In  the  form,  number  of  failures  (minimum  of  one)  and  total  operating 
time.  Some  computational  methods  to  facilitate  the  determination  of 
the  lower  confidence  bound  at  a  given  level  are  provided  In  the 
important  case  of  a  series  system.  A  numerical  comparison  Is  made 
between  the  nominal  confidence  level  and  the  actual  confidence  level, 
by  counting  the  number  of  times  the  true  reliability  la  caught  within 
the  confidence  bound  determined  by  the  likelihood  ratio  method  using 
repeated  computations  from  random  numbers  generated  with  a  specified 


distribution. 


The  problem  of  establishing  confidence  llmlte  for  the  reliability 
of  systems  has  now  extended  over  a  decade.  The  first  results  were  those 
of  Buehler  |2|  In  1957.  The  problem  he  considered  was  equivalent  with 
finding  exact  confidence  limits  for  two  components  In  aeries  with 
binomial  da 'a  on  each  component.  The  construction  of  tables  of  exact 
bounds  for  up  to  three  component  series  systems  with  various  sample  sixes 
for  tl.e  components  was  done  by  Llpow  and  Riley  (7).  This  work  In  two 
volumes  was  published  by  the  Defense  Documentation  Center. 

However,  the  bulk  of  the  tables  for  even  such  small  numbers  of 
components  made  the  use  of  simpler  approximate  confidence  limits  quite 
appealing.  Madansky  in  (9]  utilized  the  asymptotic  distribution  of  the 
likelihood  ratio  and  the  usual  practice  of  inverting  a  test  to  obtain  a 
confidence  bound,  to  yield  approximate  confidence  bounds  for  series, 
parallel  and  series-parallel  systems. 

Recently  Lentner  and  Buehler  (6)  used  the  Lehmann-Scheffd  theory 
of  exponential  families  to  find  exact  confidence  limits  for  the  specific 
case  of  components  in  series.  Due  to  the  difficulty  in  computing  these 
limits,  El  Mawaziny  and  Buehler  (4)  give  an  approximation  to  this  exact 
solution  for  the  case  where  the  sample  sizes  for  all  components  are 
large  and  the  failure  law  for  each  component  is  exponential. 

El  Mawaziny  and  Buehler  make  no  numerical  comparisons  of  their 
approximate  confidence  intervals  with  those  obtained  by  the  exact  method, 
or  the  approximations  using  the  asymptotic  distribution  of  the  maximum 
likelihood  estimates  or  the  likelihood  ratio.  They  do  show  that  under 
certain  conditions  on  the  component  samples  sizes  their  confidence 
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llmlts  are  oaymptot lcully  equal  Co  both  Cha  maximum  llkolihood  limit* 
and  to  tiie  apodal  caae  of  Bayealan  limit*  where  the  prior  "denaity"  la 
a  particular  nonprobnblllatlc  form.  They  atate  that  they  do  not  b*ll*v* 
that  their  method  Is  applicable  to  other  than  eerlea  ayeteme. 

Approximate  confidence  interval*  for  the  reliability  of  any  ayatem 
(or  atructure)  which  can  be  repreaented  by  a  monotone  Boolean  function 
of  Bernoulli  variate*  were  obtained  by  Myhre  and  Saundera  [10].  In  [10] 
the  component  failure  data  were  the  outcome  of  a  number  of  Bernoulli 
trlala  for  performance  or  nonperformance.  That  paper  waa  an  extenalon 
of  the  results  of  Madanaky  loc.  cit.»  for  aerie*  eyatema  in  that  it 
depended  upon  the  adequacy  of  the  aaymptotlc  distribution  of  the 
likelihood  ratio.  Here  we  will  follow  the  same  general  lines  of 
argument  used  in  1 10 ]  to  obtain  approximate  confidence  intervals  for 
system  reliability  from  samples  cf  component  life  lengths  but  the 
assumption  of  binomial  data  on  performance  of  each  component  la  replaced 
by  the  assumption  of  component  life  length  being  exponentially  distributed. 
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1 .  The  Cenvral  Coherent  System 

Let  the  number  of  components  in  a  given  system  be  m. 

The  state  of  the  components,  at  any  given  time  t  *  0,  Is  the 
random  vector  Jf(t)  *  (Yj (t), . . , ,YB<t))  where  Yj(t),  a 
Bernoulli  random  variable,  is  the  Indicator  of  performance  for 
the  Jth  component  at  that  time.  It  is  assumed  that  the  system 
has  a  unique  represehtatlon  as  a  nondecreasing  Boolean 
function  :,  the  functional  value  of  which,  v (£(*))•  the 
Indicator  of  the  state  of  the  system.  We  assume  without  loss  of 
generality  that  each  component  of  4  is  essential,  p.  64,  (1). 

The  reliability  of  the  jth  component  is  EY^(t)  and 
similarly  the  reliability  of  the  system  is  E1’(Y(t))  at  any 
time  t  •  0. 


If  the  life  length  of  the  component 
is  exponentially  distributed  for  j»l,...,m, 
is 


Xj  -  sup{t  >  0:  Yj (t)  ■  1} 
then  the  density  of 


(1.1)  fj(t)  ■  Aj  exp(-A^t)  for  t  >  0. 


Let  A  -  (A^,...,A  )  be  a  point  in  the  parameter  space 


CH 


{(A,,,.,  A  ) :  0 
1  m 


<  A  , 


j-1, . . . ,m} . 
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Sliue  wt-  .ire  primarily  Interested  in  Che  aystem  reliability 
at  a  proHrr Ihod  time  (in  certain  caiei  called  Che  miaalon  length) 
we  inn  w i tiii’u t  loss  of  generality  take  it  to  be  unity.  We  now 
wet 

M)  -  K !  |  Y  ( 1 )  1 

m  1-y. 

-  1  :<y)  M  exp(-»  y.)|l  -  exp(-v  ))  J 
y  J-l  J  1  J 


where  y  -  (y y  )  is  a  vertex  of  the  m-dimensional  unit 
i  m 

liypercube  and  the  wummation  ia  over  all  such  vertices.  One  sees 


(1.2)  hO)  -  N  (y)  fr  exp(->  )  +  (1-y  )  [  1  -  exp  (->  . )  ] } . 

y  J  J  J  J 

Suppose  that  Oj  _  1  identical  replications  of  the 
component  are  tested  for  j*l,...,m. 

Consequently,  observations  are  made  on  independent  random 
variables  identically  distributed  as  Xj ,  which  by  our  convention 
are  life  lengths  expressed  as  multiples  (possibly  less  than  one) 
of  the  given  mission  length.  Call  these  variables  for 

i*l , . . . i n j « 

We  postulate  that  often  in  practice  the  random  variable  which 
is  actually  observed  is  minfX^j.Z^J  where  is  the  random  time 

(again  expressed  in  terms  of  mission  length)  at  which  the  test  of  the 
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ll1'  replication  of  the  j6*’  component  will  terminate  for  any  reason 
oilier  than  failure  of  the  component.  Perhaps  Z  ^  is  degenerate  at  the 
fixed  time  at  which  the  test  is  terminated,  this  time  being  determined  by 
t lie  funds  available  for  testing  and  known  prior  to  the  test's  commencement, 


Alternatively,  could  be  the  random  time  of  failure  of  ancillary  teBt 

equipment  which  causes  the  abortion  of  the  test,  or  any  one  of  a  hundred 
other  causes.  The  important  point  is  that  we  observe  the  events 


(1.3) 


lXiJ  "  xij  1  °  lXij  <  Zij1  1-1 . 8J 


l zi j  *  Zljl  0  lXij  I  Zij 1  i"8j+l,  • • • ,Hj 


where  (the  number  of  complete  life  observations  of  the  J 

component)  is  an  observed  value  of  a  random  variable  determined  by 
the  joint  distribution  of  Z ^  i«l,...,nj,  which  is  unspecified. 

(Under  special  circumstances  the  distribution  of  the  number  of  observed 
failures  is  known  to  be  Poisson.) 

Sometimes,  of  course,  it  is  decided  prior  to  the  start  of  the 
test  of  the  jth  component  that  it  will  be  terminated  at  the  8jt^  failure 
Here  s^  would  be  a  known  parameter  and  the  time  of  completion  of  the 
test  would  become  a  random  variable.  The  required  joint  distribution 
of  Z^  i»l,...,nj  can  be  surmised  for  this  case. 

We  will  now  state  a  well-known  result  as 

Lemma  1:  If  (Z^,...,Zn)  is  a  vector  of  non-negative  random 
variables  independent  of  (X^,...,Xn)  which  are  themselves  identically 
distributed  non-negative  random  variables  with  conmon  density  f,  then 


T±=2 


t lie  likelihood  of  the  event 
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k  n 

d.4)  n  |x1_zlnxi  -xt]  n  [x.  >  zi][z1  -  *4l 

1*1  i-k+1 


is  of  the  form 

k  n 

k  TT  (<*<)  TT  (i  - 

i-1  1  k+1  3 


where  K  does  not  depend  upon  f. 

Proof.  Let  C,  be  the  joint  distribution  of  ),  then 

the  probability  of  the  event  specified  in  (1.4)  is 

k  n 

IT  f<x.i)  TT  U  "  F(z  )  ]dG(z  . z  ) 

i-1  1  i-k+1  l  i  n 

which  upon  simplification  shows  that 


Zj  'Xj  i-1, . . . ,k} 


/ 


K 


>X 


i 


JdG (z. ,  • .  • 

,k  r 


» 


z 

n 


), 


and  thus  for  any  observed  failure  times  (x^, . . . ,x^)  ,  K  does  not 
depend  upon  f .  | 

This  lemma,  in  equivalent  forms,  was  given  earlier  by  Herd  [5]  and 
Sampford  [11].  Also  in  this  connection  see  the  discussion  by  Cohen  in 
[2]. 


Throughout  what  follows  we  shall  assume  that 
have  a  distribution  such  that  with  probability  one 


for  i«l,...,n 


r 


(1.5) 


0<  Sj  IV 


It  is  clear  now  that  the  log-likelihood  function  for  the  ».u 


from  the  component  may  be  written  as 


»xp{L*(A)}  *  K.  ITf.fe,,)  TT  u-Mz„)) 
3  3  i«i  J  i-Sj+l  J 


where  is  a  constant  not  depending  upon  the  parameter  A  of  ih*j 

distribution. 

r 

From  the  specific  assumption  made  in  (1.1),  by  setting  L*  -  l*, 

j*'  1 

we  find 


(1.5.1) 


L* 


m 

(A)  -  £  (in  K.  +  a  in  A.  -  X  t  ] 

J  J  J  J  J 


where  for  typographical  convenience  we  denote  the  total  test  time 
(in  fractions  of  the  mission  length)  of  the  jth  component  by 


s 


(1.6) 


n 


J 

i-1  ij  i-Sj+l  1J 


The  logarithm  of  the  likelihood  ratio,  say  L(r),  is  given  by 


L(r)  ■  sup  L*  (A)  -  sup  L*(A). 

{ A :h (A)-r )  A  ejf 

We  now  follow  the  usual  method  of  inverting  a  test,  in  this  case 
the  likelihood  ratio  test,  in  order  to  obtain  a  confidence  interval. 

Proceeding  as  we  have  done  in  [10]  we  utilize  Wilks'  theorem  [14] 
on  the  asymptotic  behavior  of  the  logarithm  of  the  likelihood  ratio  to 
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obtain  a  confidence  set  of  level  y  for  the  system  reliability  at 
the  mission  length.  This  is 


(1.6.1) 


{r:  -2L (r)  <  x  (D> 
_  y 


2  th 

where  .  (1)  is  tlu*  quantile  of  the  Chi-square  distribution  with 

i 

one  degree  of  freedom. 


Since  the  maximum  likelihood  estimate  of  Aj  is 


'  X  m  t  J“1 , . . . ,m 

J  J 


we  see 


!.*(•)  3  sup  !.*(*) 


^  |«n  K  +  s  *n  s  -  s  (JLn  t.  +  1)]. 

,  _  1  J  .1  J  J  J 


To  maximize  !.*(•)  subject  to  the  restriction  h(A)  ■  r  we  use 
a  Lagrange  multipler  ‘,  take  partial  derivatives,  and  equate  to  zero. 


—  -  l L* ( * )  -  ih < A ) )  -  0. 


This  is  equivalent  with 


m 

1  !«,  ‘  n  f  -  ‘  t  -  5h  (A ) )  -  0, 

|-i  J  J  J  J 


ami  thus  we  obtain 
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(1.7) 


!l. 

XJ 


63jh  (A) 


where  S^h  is  the  partial  derivative  of  h  with  respect  to  A j . 

The  existence  of  3^h  follows  from  the  definition  of  h.  For  given 
6  denote  the  vector  solution  of  (1.7)  by  A (6),  assuming  presently 
that  it  exists  and  is  unique  within  which  we  shall  later  prove 
in  Theorem  3.  Note  that  A(0)  ■  A. 

Since  Lagrange  multipliers  are  being  used  it  will  be  shown  that 
the  confidence  set  may  conveniently  be  written  in  terms  of  the 
multiplier  5  rather  than  in  terms  of  the  reliability  r. 

First  define 


(1.8)  A  (6)  -  L* [ A  (5 )  ]  -  L*(A) 

for  those  values  of  6  for  which  A  (6)  exists.  It  will  be  shown  later 
in  Theorem  4  that  the  existence  of  unique  Aj(6)  implies  the  existence  of 
the  derivative  Aj(6).  Thus  it  is  possible  to  differentiate  A  with 
respect  to  6.  Since 

m 

A (<5 )  -  £  [s  in  A  (6)  -  A  (6)t  -  s.  in  s  +  s  in  t  +  s  ] 

JJ  JJJjJJJ 

m  8 

A'<5)"£ 

which  by  equation  (1.7)  shows 


A '  (5) 


53jh[A(6)]Aj(6)  -  6  ~  h[A(6)] 
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I 

From  this  equation  there  follows 

Lemma  2:  At  each  a  for  which  A ( 6 )  exists  h A  ( 6 )  is  monotone 
decreasing  in  A  if  and  only  if  A (6 )  is  monotone  decreasing  in  5 
for  ;  •  0  and  monotone  increasing  in  6  for  &  <  0. 

Using  this  result  we  obtain 

Theorem  1:  If  h'(*)  is  monotone  decreasing  across  an  interval  j 

['  ,5  ]  where  6  0  <  £+  are  two  values  of  6  for  whicn 

A  -  -4  x2d) 

Y 

then 

(1.8.1)  i  r :  h>  (-* ”)  r  •  h\  (6+)}  -  {r:  2L(r)  <_  x*(l)}  • 

l'roof.  Since  Ii-(a)  is  a  one  to  one  transformation  for  each 

2 

given  r  for  which  2L(r)  <_  there  exists  a  unique  6r  such  that 

h  \  ‘  r )  =  r.  Thus  by  definition  we  have 

L(r )  -  L*[A(6r)]  -  L*  (A  )  -  A  («r> 

and  the  parameterization  of  L  by  r  may  be  replaced  by  that  of  A  by 

To  complete  this  argument  we  must  show  that  there  exists  values  of  6 

for  which  1 ( ‘ )  exists  uniquely,  and  hopefully  can  be  easily  found, 

and  also  that  there  exists  values  of  6  in  an  interval  about  zero  for 

which  h » ( • )  is  decreasing.  We  turn  to  the  first  task  now. 


i 
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Tht>  transformation  A  from  31  into  for  fixed  6  we  define 


.th 


by  setting  its  j  component  for  j-l,...,m 


(1.9) 


Aj(X;4)  "  tj+sa^h (A)  * 


This  is  suggested  by  solving  (1.7)  as  if  9 ^ h  were  constant. 

T.t  is  clear  that  Aj(A;6)  is  continuous  in  6.  Since  9jh(A)  <  0, 
see  equation  (1.14)  below,  A^(A;*)  is  an  increasing  function  of  6 
for  tj  >  -6D jh (A ) . 

We  now  quote  from  [10]  a  result  easily  proved,  that  for  a,b  e  (0,1) 
and  <j,B  any  real  numbers 


(1.10) 


aa-b8  |  <_  |  i-b  |  +  |  a-3 1 . 


For  a  given  structure  $  define  the  criticality  of  the  j 
component  in  the  structure  by 


th 


(1.11) 


cj  *  X  *<y|J:1)  "  ®(y|j:0) 


where  (yjjsx)  -  (y^ » • • • »yj_^ ,x,yj+T * * '  ’ » y^  ^  for 
We  now  have 


x  *  0,1. 


Theorem  2:  For  each  3>  (or  h)  and  all  6  such  that 


m  sc  . 

(1.12)  min(t.  , . . .  ,t  )  >  6  >  0  and  7  — ^  ■*  ■■  <  -f 

1  "  1-1  (tj-6)2  6 


(1.13) 


m  sc  . 

•  0  and  ^  -M"  <■  7-TT  , 

w  tj2  I5' 


the  transformation  A(  is  a  contractive  map  of  the  complete  metric 

space  C 7f,d )  into  Itself  where 


1 

d('  ,i  )  =  -  ^  |  |  for  \,p 


Proof.  It  must  be  shown  that  there  exists  p  c  [0,1)  such  that 


d [  A ( •  )  ,A(n  )  ]  pd  (A  ,u)  for  all  A  ,p  c  3i. 


Note  that  throughout  this  proof  we  shall  omit  showing  the  dependence  of 
both  A  and  p  upon  ■'  in  order  to  simplify  the  notation.  Now 


s  .  5  [  D  .h  (p )-9  h  (a  )  ] 


AJ(')  Ai(u)  =  [tj+{31h(A))[tj+5^jh(p))  * 


Since  by  known  properties  of  coherent  systems,  see  [1] 


-Hy)  =  y  j  *'  (y  I  j :  i )  +  ( i-y  j  )(f  (y  I J  :0> 


substitute  Y^(l)  for  and  take  expectations  to  obtain 


-A  -A 

h(»)  =  e  Jh(A (j : 1)  +  (1-e  J)h(A  | j  :0) . 


:-.h(\)  =  e  J[h(A|j:0)  -  h(A|j:l>). 


(1.14) 
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If  we  continue  lie  expansion  we  obtain 

f 

h(A)  =  e  Y  {(Ky|j:0)  -  4>(y|j:l)]  [y.e  1  +  (l-y.)(l-e  i)]}. 

1  y“  W  1  1 

Tlius  writing  out  [jjh(v)  -  3jh(\);  using  the  definition  of  (1.11)  and 
applying  the  inequality  (1.10)  repeatedly  on  the  last  fraction  we  find: 

™  -I  4 

I  a  h  (p )  -  3  h(A)|  <_  c  2  |e  -  e  1|. 

J  J  J 

^  i  ~ ^  i 

Since  |e  -  e  |  <^  | 

|  3jh  (p )  -  3  h(A)|  <_  mCjd  (A  ,p ) 

and 

m  s  c  5d(A,y) 

d  (A(A  )  , A(u ) )  <  ^  [tj+fiajMAJHtj+fidjh^)]  ‘ 

Since  each  component  is  essential  it  follows  from  (1.14)  that 

3  h(A)  <  0  as  0  «'  h(A  |  j  :1)  -  h (X  )  j  :0)  <  1.  It  also  follows  from 

"A1 

(1.14)  that  Djh(A)  >_  -e  J  >_  -1.  Thus  it  is  sufficient  to  require 
that  (1.12)  hold  for  6  >  0  and  (1.13)  hold  for  6  <  0.|| 

For  any  A 0  c  Jf  we  define  the  sequence 

An(6)  -  A(An_1(6)  ;6)  n-1,2,. .. 

where 


A°(6)  -  A0. 
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Theorem  3:  For  every  6  in  the  neighborhood  of  zero  defined  by 
the  inequalities  (1.12)  and  (1.13),  a  unique  solution  to  the  system  of 
equations  (1.7)  exists,  call  it  A  (6).  It  can  be  found  for  any  initial 
point  A0  r  31  as  lim  A°(6)  ■  A(6). 

n  ►. 

Proof.  That  A ( • ; 6 )  for  6  within  the  prescribed  neighborhood 
of  zero  has  a  unique  fixed  point 

A  (6 )  =  lim  An(6)  -  A(A  (6 ) : 6 ) 

follows  from  the  known  behavior  of  contractive  maps,  e.g.  see  p.  27, 

[8].|| 

Thus  we  have  established  the  first  claim  which  was  made,  namely 
that  A (6)  exists  uniquely.  Now  we  prove 

Theorem  4:  The  function  A(*)  is  a  continuously  differentiable 
function  within  the  neighborhood  of  zero  prescribed  by  (1.12)  and  (1.13). 

Proof.  Fix  5  within  the  prescribed  neighborhood  and  let  B  be 
the  vector  valued  function  with  its  j1”*1  coordinate  defined  by 

B  (A  : 6 )  =  69  ,h(A  )  +  t  -  r1  . 
j  J  J  j 

By  Theorem  3  the  equation  B(A:5)  *  0  has  a  unique  solution,  call  it  A. 
Thus  the  Jacobian  is  not  zero,  i.e. 

det  9 ±Bj (A ; 6 )  +  0. 

From  the  implicit  function  theorem  the  continuous  differentiability  of 
the  function  B  is  inherited  by  A.  Thus  if  A  exists  it  is  continuously 
differentiable. | | 
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In  view  of  Theorems  3  and  4  it  remains  only  to  show  that  hA(6) 
is  a  decreasing  function  of  6.  We  then  have  the  result  that  for  any 
coherent  structure  an  approximate  confidence  interval  for  the 
reliability  function  !i  is  given  by  (1.8.1)  in  Theorem  1. 

Theorem  5:  For  any  coherent  structure  $  with  reliability 
function  h  and  any  failure  data  such  that  s  1,  there  exists 
a  neighborhood  of  zero  in  A.  across  which  hA(6)  is  decreasing  in  6 


Proof.  Take  the  derivative  of  both  sides  of  (1.7)  with  respect 
to  6,  primes  denoting  such  differentiation,  we  obtain  for  j*l,...,m 


(1.15) 


')  hA  +  6  V  3  h  (A  )  •  A  ' 

J  i-1  iJ  1 


where  we  have  omitted  the  argument  A.  Multiply  by  Aj  on  both  sides 
of  (1.15)  and  sum  to  obtain 


V  _ L_  pi%2 

^  O  V  1  ' 

3-1  <y 


1  \h(A) 

j  =  l  J 


m  m 

i;  +  5  i  i  3 

3  i-1  j-1  J  J 


which  is  equivalent  with, 


(1.16)  “  d5  h[A ( 


m  s 

'01  =  ^  (a!)Z  +  6  £  £  3 

j-i  dj)2  J  i-i  A 


± jh  (A)A^A j . 


For  6  =  0,  (1.16)  is  positive  and  from  the  continuity  of  77  h(  A  (6)1 

do 

it  follows  that  there  exists  a  neighborhood  of  6  about  zero  such  that 
h\(A)  is  a  decreasing  function  of  5. | J 


Thus  for  any  coherent  structure  where  the  data  on  each  component's 
reliability  consists  of  knowing  the  times  of  failure  or  the  lengths  of 
time  a  component  worked  without  failure,  assuming  the  life  length  of 
the  individual  components  are  exponentially  distributed,  it  is  possible 
to  find  an  approximate  confidence  interval  of  some  level  about  the 
..laximum  likelihood  estimate  of  the  system  reliability  at  the  given 
(mission)  time.  Of  course,  it  may  not  be  possible  that  an  approximate 
confidence  interval  of  arbitrary  level  can  be  found  for  every  structure 
and  every  set  of  sample  da. a,  but  it  is  always  possible  to  obtain  an 


interval  estimate. 
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2.  The  Series  System 

For  any  coherent  structure  with  associated  reliability  function  h, 
there  exists  a  series  system,  the  reliability  of  which  bounds  the  system 
reliability  h  from  below,  i.e. 


(2.1) 


h(A)  >_  exp 


m 

-  -  Vi 
J-l  2  2 


where  Vj  is  the  number  of  minimal  cuts  of  h  (or  4>)  which  contain  j 
as  a  member.  For  the  exact  definition  of  minimal  cut  and  a  more  detailed 
discussion  of  this  point,  see  [12]. 


Because  of  this  fact  we  will,  in  this  section  consider  only  the  case 
£n  h(A)  ■  “EvjAj  and  obtain  a  lower  confidence  bound.  We  think  this  it. 
the  most  important  case. 


We  will  follow  the  general  procedure  which  has  been  presented  in 
the  preceding  section.  However,  as  a  convenience  we  reparametrize  by 
Imposing  the  restriction  in  the  form  that  the  £n  h(A)  is  constant, 
and  equate 

[L*(A)  -  <5  £n  h  (A )  ]  -  0. 


Thus  we  obtain,  instead  of  (1.7),  the  equation 


(2.2) 


-5uj 


which  has  the  explicit  solution 


a. (6)  -  T~t~ 

j  tj-SVj 


• • • ,m. 
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Thus  in  this  particular  case,  it  is  unnecessary  to  use  the  general 
contractive  map  theorems  in  order  to  calculate  x(fi)  ■  (j^(6)» • • • » Xm ( 6 ) ) • 
For  this  special  case  we  find,  using  (1.5.1)  and  (1.8), 

"l  T  .  T, 

A (6)  -  1  s  [1  -  — +  Jln( — “)] 
j-i  3  tj  5  Tj  { 

where  tj  ■  t^/v^  f°r  j“l,...,m.  In  accord  with  finding  only  a  lower 
confidence  bound,  we  want  lo  solve  for  x  >  0  such  that 


(2.3) 


A (x)  -  -h  x2d>. 

Y 


where  is  the  100yth 
with  one  degree  of  freedom. 
Then 


percentile  of  the  Chi-square  distribution 

Call  the  solution  of  (2.3)  the  value  6+. 

Y 


h( A (6  ) ]  -  exp 


m 

-  I 


J-l  ^ 


-6 


is  an  (approximate)  lower  bound  for  h(A)  of  level  ,  rather  than 

of  level  y i  since  we  are  obtaining  only  a  one-sided  confidence  bound. 


We  now  exhibit  a  practical  method  for  the  determination  of  6+. 

Y 

For  a  given  y,  0  <  y  <  1  there  exists  an  x  >  0  such  that  f(x)  ■  0 
where 


f  (x) 


-A(x)  -  h  x*(l) 


m 

2 


Sj  -  H  X*<1). 


(2.4) 


The  solution  of  f(x)  *0  must  be  unique  since 


m  s 

f'(x)  -  X  y  - i-r  >  0 

i-i  (Tt-xr 


whenever  0  <  x  <  min(T^, . . .  ,-r^)  ■  which  is  a  more  stringent 

condition  than  (1.12).  But  note  that 


m  t  .+x 

f"(x)  -  Is.  — - r  >  0 

i-1  (xi-x)J 

for  0  <  x  <  t 

Thus  we  see  that  both  f',f",  which  are  continuous,  do  not  vanish 
for  0  <  x  <•  t^.  These  are  sufficient  conditions  that  the  Newton 
iteration  procedure,  namely 


f  (x  - ) 
n-1 

x  ■  x  ,  -  ttt - r 

n  n-1  f'(x  ,) 
n-1 


n-1, 2 


, . . . , 


will  converge  to  the  value 


We  now  state  the 


Theorem  6:  Let  the  data  i-l,...,m  be  given,  where  t^ 

represents  the  total  test  time  for  the  ith  component  and  s^  >_  1, 
represents  the  number  of  failures  of  that  component  during  that  time. 
If  v ^  is  the  criticality  of  the  ith  component  in  any  structure  being 
measured  by  the  number  of  minimal  cuts  of  the  system  which  contain  the 
i^  component  as  a  member,  and  x  -  t^/v^,  then 
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(2.5) 


exp 


_  V  _ 1 _ 

J-l  (Tj/t )-&+ 


for 


0  <  t  < 


T(l) 


0 


for  t  > 

-  (1)  Y 


Is  a  lower  confidence  bound  of  level  ,  where  0  <  y  <  1  is 

preselected,  for  the  system  reliability  for  all  time  t  >  0.  The  quantity 

i  is  the  solution  of  the  equation  f (x)  ■  0  and  f  was  defined  in 

*1 

equation  (2.4). 

Proof.  We  have  only  to  remark  that  the  selection  of  a  time  scale 
in  terms  of  mission  length  was  arbifary  so  that  the  bound  derived 
previously  can  be  used  for  the  reliability  at  all  times.  |  I 
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3.  A  Numerical  Illustration 

As  a  numerical  Illustration  of  the  behavior  of  the  previously 
outlined  procedure  we  shall  study  its  performance  on  a  specified  series 
system  with  ten  components.  Take  m  ■  10.  We  shall  assume  that 
Aj  •  1  for  j«l,...,10.  Thus  the  true  reliability  of  the  system  is, 
by  (1.2)  e  for  t  >  0.  We  also  assume  that  nj  “  vj  “  ^  anc* 

thus  t.  ■  t.  for  j*l,...,m. 

J  J  J 

We  will  now  generate  10  exponential  variates,  t^ . t  ,  each 

with  unit  mean  and  solve  the  equation  f(x)  ■  0  where 


(3.1) 


f  (x) 


C 

Y 


and  is  a  constant  determined  by  the  nominal  level  of  the  confidence 

bound . 

We  find  that  for  y  -  .95,  C  9$  -  11.353.  Call  the  solution  of 

(3.1)  the  value  6*.  The  lower  bound  for  the  reliability  is  by  (2.5) 

m  . 

(3.2)  exp  -  X  ( - iu“T)  .  for  0  <  t  <  6+t,.., 

i-1  t . / t-6+  Y  (1) 

i  Y 

and  the  confidence  level  is  97.5. 

Generating  forty  independent  observations  of  an  exponential  variate 
with  unit  mean,  resulted  in  the  following  four  samples: 
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(1) 

(2) 

(3) 

(4) 

1.305 

1.889 

.2331 

3.258 

1.346 

3.864 

.8881 

.1848 

1.483 

.7789 

.04059 

.7594 

.2931  .4825  .8329  2.420 

.2282  1.171  .1948  .8638 


.3489 

.5642 

1.332 

.5307 

1.342 

.3525 

2.105 

.06695 

.7171 

.9037 

1.144 

.8239 

.3397 

.5590 

1.203 

.2452 

1.126 

.5445 

.07427 

.2306 

Solving  (3.1)  by  machine  program  and  calculating  (3.2)  for  each 
set  of  observations  yielded  the  confidence  bounds  which  are  summarized 
visually  in  Figure  L. 


Figure  1.  97.5%  lower  confidence  bounds  on  the  reliability  of  a  10 

component  system  based  on  independent  samples. 

(0)  is  the  true  reliability  e  for  t  >  0 

(i)  is  the  confidence  bound  based  on  sample  (i) 
for  i=l,2,3,4. 
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As  a  further  illustration  of  the  stochastic  behavior  of  this 
method  we  shall  repeat  the  entire  procedure  two  thousand  times  but 
instead  of  calculating  the  entire  lower  bound  function  we  shall  only 
compute  (3.2)  at  t  ■  .01.  Then  we  shall  make  a  frequency  histogram 
of  the  values  of  the  reliability  bounds  for  the  specified  value 
t  ■  .01.  The  twenty  percentile  points,  successive  differences  being 
.05,  of  the  empiric  distribution  are: 


Percentile 

Percentile  Point 

.05 

0.106621 

.10 

0.297174 

.15 

0.417770 

.20 

0.492659 

.25 

0.552798 

.30 

0.596595 

.35 

0.632835 

.40 

0.665383 

.45 

0.689521 

.50 

0.707729 

.55 

0.730936 

.60 

0.750021 

.65 

0.767430 

.70 

0.782541 

.75 

0.798792 

.80 

0.815686 

.85 

0.833291 

.90 

0.851064 

.95 

0.877483 

1.00 

0.940218 

A  graph  of  the  empiric  distribution  is  given  in  Figure  2.  Note 
that  the  true  reliability  of  e  -  .904837  slightly  exceeded  the 
nominal  97.5  percentile.  In  fact  the  actual  count  was  1979  values 
less  than  e  out  of  the  2000  observations  sampled.  This  was 


approximately  the  99th  percentile. 


Cumulative  percentage  less  than  or  equal 
to  the  reliability 


Reliability 


Figure  2. 


Empiric  distribution  of  the  97.5  per  cent  lower  confidence 
bound  using  2000  observations  when  the  true  reliability  is 
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Thus  we  see  the  stated  confidence  level  of  97.5  percent  appears 
to  be,  in  this  instance,  slightly  conservative.  To  check  this  we 
repeated  this  entire  experiment  a  second  time  and  we  observed  1990  out 
of  2000  observations  of  the  reliability  bound  were  less  than  e  ' ^ 
which  increases  our  suspicion  of  a  slightly  conservative  tendency 
for  the  level  of  confidence  in  this  case. 


-27- 


References 

[1]  Birnbaum,  Z.  W.,  Esary,  J.  D. ,  Saunders,  S.  C.  ,  Multi-component 
Systems  and  Structures  and  Their  Reliability,  Technometrics 
55-77,  1961. 

[2]  Buehler,  R.  J.,  Confidence  Intervals  for  the  Product  of  Two 
Binomial  Parameters,  J.  Amer.  Statist.  Assoc.  5J2,  482-493,  1957. 

[3]  Cohen,  A.  Clifford,  Life  Testing  and  Early  Failure,  Technor.e  tries 
539-545,  1966. 

[4]  El  Mawaziny,  A.  H.  and  Buehler,  R.  J.,  Confidence  Limits  for  the 
Reliability  of  Series  Systems,  J.  Amer.  Statist.  Assoc.  62 _,  1452-1459, 
1967. 

[5]  Herd,  G.  R. ,  Estimation  of  the  Parameters  of  a  Population  from  a 
Multi-censored  Sample,  Ph.  D.  Dissertation,  Iowa  State  College,  1956. 

[6]  Lentner,  M.  M.  and  Buehler,  R.  J.,  Some  Inferences  about  Gamma 
Parameters  with  an  Application  to  a  Reliability  Problem,  J.  Amer. 
Statist.  Assoc.  5j8,  670-677,  1963. 

[7]  Lipow,  M.  and  Riley,  J.,  Tables  of  Upper  Confidence  Limits  on 
Probability  of  Failure  for  1,  2,  and  3  Component  Serial  System 
(in  two  volumes),  Space  Technology  Laboratories  TR-59-0000-00756, 

1959. 

[8]  Liustemick,  L.,  Sobolev,  V.,  Elements  of  Functional  Analysis, 

Ungar  Publishing  Co.,  1961. 

[9]  Madansky,  Albert,  Approximate  Confidence  Limits  for  the  Reliability 
of  Series  and  Parallel  Systems,  Technometrics  4.,  495-503,  1965. 


-28- 


[10]  Myhre,  J.  M.,  Saunders,  S.  C.,  On  Confidence  Limits  for  the 
Reliability  of  Systems,  Am.  Math .  Statist.  22.,  1968. 

[11]  Sampford,  M.  R. ,  The  Estimation  of  Response-Time  Distributions, 

Part  II,  Biometrics  9.,  307-369,  1952. 

[12]  Saunders,  S.  C.,  On  Confidence  Limits  for  the  Performance  of  a 
System  when  Few  Failures  are  Encountered,  Boeing  Scientific  Research 
Laboratories  document  Dl-82-0676,  November  1967;  submitted  for 
publication  in  J.  Amer.  Statist.  Assoc. 

[13]  Springer,  M.  D.  and  Thompson,  W.  E. ,  Bayesian  Confidence  Limits 
for  Reliability  of  Redundant  Systems  when  Tests  are  Terminated  at 
First  Failure,  Technometrics  10,  29-36,  1968. 

[14]  Wilks,  S.  S.,  The  Large-sample  Distribution  of  the  Likelihood 
Ratio  for  Testing  Composite  Hypothesis,  Ann.  Math.  Statist.  _9, 

60-62,  1938. 


